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Abstract 

The synthesis of the negative-strand [(-)-strand] complement of the —30 kilobase, positive-strand [(+)-strand] coronaviral 
genome is a necessary early step for genome replication. The identification of c/s-acting elements required for (— )-strand 
RNA synthesis in coronaviruses, however, has been hampered due to insufficiencies in the techniques used to detect the (— 
)-strand RNA species. Here, we employed a method of head-to-tail ligation and real-time quantitative reverse transcription 
polymerase chain reaction (qRT-PCR) to detect and quantitate the synthesis of bovine coronavirus (BCoV) defective 
interfering (Dl) RNA (— ) strands. Furthermore, using the aforementioned techniques along with Northern blot assay, we 
specifically defined the c/s-acting RNA elements within the 3'-terminal 55 nucleotides (nts) which function in the synthesis 
of (— )- or (+)-strand BCoV Dl RNA. The major findings are as follows: (i) nts from -5 to -39 within the 3'-terminal 55 nts are 
the c/s-acting elements responsible for (— )-strand BCoV Dl RNA synthesis, (ii) nts from —3 to —34 within the 3'-terminal 
55 nts are c/s-acting elements required for (+)-strand BCoV Dl RNA synthesis, and (iii) the nucleotide species at the 3'-most 
position (—1) is important, but not critical, for both (— )- and (+)-strand BCoV Dl RNA synthesis. These results demonstrate 
that the 3'-terminal 55 nts in BCoV Dl RNA harbor c/s-acting RNA elements required for both (-)- and (+)-strand Dl RNA 
synthesis and extend our knowledge on the mechanisms of coronavirus replication. The method of head-to-tail ligation and 
qRT-PCR employed in the study may also be applied to identify other c/s-acting elements required for (— )-strand RNA 
synthesis in coronaviruses. 



@ PLOS | one 



Citation: Liao W-Y, Ke T-Y, Wu H-Y (2014) The 3'-Terminal 55 Nucleotides of Bovine Coronavirus Defective Interfering RNA Harbor Gs-Acting Elements Required 
for Both Negative- and Positive-Strand RNA Synthesis. PLoS ONE 9(5): e98422. doi:10.1371/joumal.pone.0098422 

Editor: Volker Thiel, Kantonal Hospital St. Gallen, Switzerland 

Received January 29, 2014; Accepted May 2, 2014; Published May 22, 2014 

Copyright: © 2014 Liao et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by grant 99-231 3-B-005-024 from the National Science Council (NSC) of Republic of China. The funders had no role in study 
design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 

* E-mail: hwu2@dragon.nchu.edu.tw 



Introduction 

After entry into a cell, the positive-strand [(+)-strand] RNA virus 
genome first serves as a template for synthesis of viral proteins 
required for subsequent virus replication. The synthesis of the 
negative-strand [(— )-strand] complement from the (+)-strand 
genome is the first step of genome replication and is assumed to 
be initiated from the 3' end of the (+)-strand RNA genome. The 3' 
end of the synthesized (— )-strand RNA genome is then employed 
as the site for the initiation of the (+)-strand RNA genome. 
Therefore, the RNA elements on the 3' end of (+)- and (— )-strand 
RNA genome may play critical roles in targeting the viral 
replication complex to initiate the synthesis of their counterparts. 

Bovine coronavirus (BCoV), a (+)-strand RNA virus in the 
subfamily Coronavirinae, the family Coronaviridae and the order 
Nidoviraks, belongs to the genus betacoranaoiruses, which includes the 
mouse hepatitis coronavirus (MHV) [1,2]. Many as-acting RNA 
elements in the coronavirus genome required for replication 
[interpreted as (+)-strand RNA synthesis] have been identified in 
the past two decades with a full-length infectious cDNA or 
coronavirus defective interfering (Dl) RNA, a surrogate for the 
coronavirus genome [3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]. 



For as-acting RNA elements in the 3' UTR of coronavirus, it 
has been demonstrated that the minimal requirement for 
replication maps 436 nucleotides (nts) from the 3' end of MHV 
Dl RNA according to deletion analyses [19,20]. Within this 
region, two higher-order structures on the upstream end of the 3' 
UTR have been identified as as-acting RNA elements that are 
essential for coronavirus replication, including a 61-nt bulged 
stem-loop in MHV [12,13] and a hairpin-type pseudoknot in 
BCoV Dl RNA [15]. Further downstream of the 3' UTR is a 
hypervariable region (HVR) which forms a complex secondary 
structure. Although this secondary structure was shown to be 
important for the efficient replication of MHV Dl RNA [14], a 
subsequent study using an intact MHV genome demonstrated that 
the secondary structure is nonessential for RNA synthesis but plays 
a pivotal role in pathogenesis [21]. Downstream of these higher- 
order structures are 3 '-terminal 55 nts and poly(A) tail. Although 
coronaviral poly(A) tail has been demonstrated to be a critical 
element for coronavirus replication [22], the role of the entire 3'- 
terminal 55 nts in the (+)-strand RNA synthesis remains to be 
determined. 
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In contrast, studies on the oV-acting RNA elements required for 
(— )-strand RNA synthesis in coronaviruses have been hampered 
by the limit amount of (— )-strand RNA present in the virus- 
infected cells and the lack of sufficient methods to detect (— )-strand 
coronaviral RNA in low numbers. For example, there are only an 
estimated ~5 molecules of BCoV DI RNA (— ) strand per cell in 
BCoV DI RNA-transfected BCoV-infected human rectal tumor 
(HRT)-18 cells at 24 hour postinfection (hpi) of virus passage 1 
(VP1) [23,24]. In order to overcome the detection deficiency and 
the false positive results of the (— )-strand DI RNA detection 
caused by artifact of (— )-strand DI RNA transcripts generated by 
T7 RNA polymerase using plasmid DNA as a template [25], a 
head-to-tail ligation method and reverse transcription polymerase 
chain reaction (RT-PCR) have been developed to detect (— )- 
strand BCoV DI RNA [26]. This method, therefore, may be 
employed to determine the requirements of ciy-acting elements for 
the synthesis of BCoV (— )-strand DI RNA. 

In this study we used head-to-tail ligation and real-time 
quantitative reverse transcription polymerase chain reaction 
(qRT-PCR) to detect and quantitate BCoV (-)-strand DI RNA 
synthesis. Using this technique along with Northern blot assay, we 
specifically defined the as-acting elements within the 3 '-terminal 
55 nts that are required for the synthesis of (+)- or (— )-strand 
BCoV DI RNA by simultaneous comparison of their efficiency of 
RNA synthesis. Our findings suggest that the 3' -terminal 55 nts 
harbor as-acting RNA elements required for both (+)- and (— )- 
strand RNA synthesis. 

Results 

The 3'-terminal 55 nts within the individual subgroups of 
betacoronaviruses were highly conserved 

Before the outbreak of severe acute respiratory syndrome 
(SARS), betacoronaviruses were considered to include one lineage. 
After the discovery of SARS-CoV, other novel subgroups of the 
betacoronaviruses have been proposed and the original group 2 
coronaviruses, which included BCoV and MHV, became desig- 
nated as 2a coronaviruses [27,28]. Many cis-acting elements in the 
3' UTR of BCoV and MHV-A59 have been characterized 
[12,13,15,22,29]. To examine whether the 3'-most 55 nts have 
conserved features within the individual subgroups of betacor- 
onaviruses, the 3'-terminal 55-nt sequence was aligned according 
to the published sequences (see GenBank Accession Nos. in Fig. 1) 
and the sequence identity within the subgroup 2a, 2b, 2c and 2d 
coronaviruses was 71%, 100%, 88% and 82%, respectively (Fig. 1). 
It is apparent that the aligned 3'-terminal 55-nt sequence was 
highly conserved (identified by shading) across the coronaviruses 
within the individual subgroups of betacoronaviruses, suggesting 
that this region may harbor ay-acting elements critical for 
coronavirus replication. 

The 3'-terminal 55 nts on BCoV DI RNA are required for 
(-)-strand RNA synthesis as analyzed by head-to-tail RNA 
ligation and qRT-PCR 

The coronavirus DI RNA has been employed as a surrogate for 
coronavirus genome to study the cis-acting elements required for 
coronavirus replication [3,5,6,11,15,20,22,29,30]. The identifica- 
tion of as-acting RNA elements that are essential for the (— )- 
strand DI RNA synthesis in coronaviruses, however, has been 
restricted by the low numbers of (— )-strand DI RNA produced 
from the DI RNA-transfected helper virus-infected cells [23,24]. 
Although RT-PCR may be used to detect the (-)-strand DI RNA 
molecules, a false positive result may occur due to the synthesis of 
the copy-back (— )-strand DI RNA transcripts generated by T7 



RNA polymerase in vitro using plasmid DNA as a template [25]. 
Thus, to overcome the detection problems caused by the low 
numbers of (— )-strand DI RNA and to avoid the detection of the 
(— )-strand DI RNA resulting from the copy-back transcripts, we, 
based on the previous study [26], (1) created a hybrid BCoV DI 
RNA construct marked with the MHV 3' UTR and 25 nts of 
poly(A) tail (Fig. 2A) to distinguish it from BCoV helper virus 
genome during RT-PCR and (2) conducted head-to-tail ligation of 
RNA and RT-PCR to differentiate the coronaviral polymerase- 
generated (— )-strand DI RNA from T7 RNA polymerase- 
generated artifact (— )-strand DI RNA (Fig. 2B), as the primers 
used in the RT-PCR cannot yield a product from the copy-back 
(— )-strand transcripts after ligation [26]. Note that the sequence in 
the 3'-termianl 55 nts of the 3'UTR between BCoV and MHV is 
highly conserved (Fig. 1) and it has been shown that the replication 
efficiency of BCoV DI RNA with the MHV 3' UTR is almost the 
same with that of BCoV DI RNA with BCoV 3'UTR [26]. 
Accordingly, to test whether the method can be applied to 
examine the requirement of the 3'-most 55 nts for (— )-strand DI 
RNA synthesis in the context of BCoV DI RNA, BM25A 
construct in which the 3 '-most 55 nts were intact and A55 mutant 
in which the 3'-most 55 nts were deleted were created (Fig. 3A). 
After transfection of BM25A and A55 transcripts into BCoV- 
infected HRT- 1 8 cells, RNA was extracted at different time points 
as indicated in Fig. 3B, head-to-tail ligated, and subjected to RT- 
PCR. An RT-PCR product with a length of -150 base pair (bp) 
was obtained from BM25A-transfected BCoV-infected cells 
(Fig. 3B, lanes 2-8) but not from A55-transfected BCoV-infected 
cells (Fig. 3B, lanes 10-16), mock-infected cells (Fig. 3B, lane 18), 
BCoV-infected cells (Fig. 3B, lane 19), BM25A-transfected mock- 
infected cells (Fig. 3B, lane 20), or a mixture of BCoV-infected 
cellular RNA extracted at 10 hpi and 200 ng of BM25A transcript 
(Fig. 3B, lane 21). These results, along with sequencing analysis 
(Fig. 3C), confirmed that the ~ 150-bp RT-PCR product obtained 
from BM25A-transfected BCoV-infected cells was specific from 
head-to-tail-ligated coronavirus polymerase-generated (-)-strand 
DI RNA. Because DI RNA may undergo recombination with 
helper virus genome under select pressure [31,32,33,34], the 
detected (-)-strand DI RNA in Fig. 3B may be derived from the 
recombinant genome rather than input (+)-strand DI RNA. To 
seek evidence of a DI RNA reporter-containing recombinant 
genome, the RNA samples extracted from the BM25A- and A55- 
transfected BCoV-infected cells at 8 hour posttransfection (hpt) 
were tested for a 1,639-nt helper virus genome-integrated 
reporter-containing sequence between the membrane (M) protein 
gene and the DI RNA 3' UTR using RT-PCR [26]. No RT-PCR 
product (Fig. 3D, lanes 2-3) was observed, suggesting that the 
detected ~ 150 bp RT-PCR products do not result from a 
recombinant helper virus genome. To compare the efficiency of 
(-)-strand DI RNA synthesis between BM25A and A55, head-to- 
tail ligated RNA extracted at 8 hpt was subjected to qRT-PCR 
with the same primer set used for the RT-PCR analysis. The 
results (Fig. 3E, left panel) suggest that the efficiency of (— )-strand 
DI RNA synthesis from A55 was significantly impaired in 
comparison with that from BM25A. Thus, we conclude that (1) 
using the head-to-tail RNA ligation and qRT-PCR allows us to 
detect and quantitate (-)-strand BCoV DI RNA and (2) the 3'- 
most 55 nts in BCoV DI RNA are required for (-)-strand DI RNA 
synthesis. 

Identification of c/s-acting RNA elements within the 3'- 
terminal 55 nts for (-)- and (+)-strand DI RNA synthesis 

We have demonstrated in the present study that the 3 '-most 
55 nts are required for the synthesis of (-)-strand BCoV DI RNA 
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Figure 1. Alignment of the 3 '-terminal 55-nt sequence within the individual subgroups of betacoronaviruses. The conserved sequence 
within individual subgroups is identified by shading. Position numbers - 1 and -55 below the sequence indicate the first and 55th nt counted from 
the poly(A) tail, respectively. Abbreviations: CoV, coronavirus; BCoV, bovine coronavirus-Mebus; OC43, human coronavirus-OC43; HECoV-4408, 
human enteric coronavirus-4408; HEV, porcine hemagglutinating encephalomyelitis virus-TN11; ECoV, equine coronavirus-NC99; MHVA59, mouse 
hepatitis virus-A59; MHVJ, mouse hepatitis virus-JHM; PV, puffinosis virus; SARS-CoV, SARS coronavirus; B-S-HKU3, bat SARS coronavirus HKU3; B-S- 
HKU4, bat SARS coronavirus HKU4; B-S-HKU5, bat SARS coronavirus HKU5; BCHKU9-1, bat coronavirus HKU9-1; BCHKU9-2, bat coronavirus HKU9-2; 
BCHKU9-3, bat coronavirus HKU9-3; BCHKU9-4, bat coronavirus HKU9-4. GenBank Accession Nos. for the sequences studied here are as follows: 
U00735 for BCoV-Mebus, AF523847 for HCoV-OC43, AF523848 for HECoV-4408, AF523849 for HEV-TN11, AF523850 for ECoV-NC99, NC001846 for 
MHV-A59, X00990 for MHV-JHM, AJ54471 8 for PV, NC00471 8 for SARS-CoV, NC009694 for B-S-HKU3, NC00901 9 for B-S-HKU4, NC009020 for B-S-HKU5, 
EF065513.1 for BCHKU9-1, EF065514.1 for BCHKU9-2, EF065515.1 for BCHKU9-3, and EF065516.1 for BCHKU9-4. 
doi:1 0.1 371 /journal.pone.0098422.g001 



using the head-to-tail RNA ligation and qRT-PCR (Fig. 3). The 
map positions of the CM-acting RNA elements required for (-)- and 
(+)-strand DI RNA synthesis within this 55-nt region, however, 
remain to be determined. To define the specific m-acting RNA 
elements within the 3'-terminal 55 nts which function in the 
synthesis of BCoV DI RNA (-) or (+) strands, a series of DI RNA 
constructs with deletions either from the 3' or 5' end of the 55-nt 
region were generated by mutagenesis (Fig. 4A) and tested by the 
methods depicted in Fig. SI and described as follows. The T7 
RNA polymerase-generated transcript of each construct was 
transfected into BCoV-infected HRT-18 cells and virus within 
the transfected cells was referred to as virus passage 0 (VPO). The 
total cellular RNA was collected at 8 hpt of RNA transcript to 
detect the synthesis of (-)-strand DI RNA by qRT-PCR. Since the 
transfected transcript was still in the BCoV-infected cells at 8 hpt 
as shown in Fig. 3E (right panel), it is difficult to differentiate 
between the coronaviral polymerase-generated (+)-strand DI RNA 
and T7 RNA polymerase-generated (+)-strand DI RNA by 
Northern blot assay at this time point. Therefore, to determine 
whether the detected (+)-strand DI RNA is specifically synthesized 
by the BCoV polymerase, the supernatant virus collected at 
48 hpt, which is now called virus passage 1 or VP1, was used to 
infect fresh HRT-18 cells (Fig. SI). The total cellular RNA 
collected at 48 hpi of VP1 was used to detect the synthesis of (+)- 
strand DI RNA by the Northern blot assay. If the T7 RNA 



polymerase-generated DI RNA, which is transfected at VPO, is not 
able to replicate, the transfected (+)-strand DI RNA cannot be 
detected at 48 hpi of VP1 [5,6,15,35,36]. The signals detected by 
Northern blot assay were quantitated and presented in percentages 
as a result of (+)-strand DI RNA synthesis [relative to the wt 
BM25A (100%)] . These results are summarized in Figs 4B and 4C. 
As shown in Fig. 4B, DI RNA constructs with a 1-, 2-, 3-, and 4-nt 
deletion (constructs Al, A2, A3 and A4, respectively) at the 3' 
terminus of the 55-nt region still maintained the synthesis of (— )- 
strand DI RNA (Fig. 4B) although the synthesis efficiency was 
impaired (—55% to —80% of the wt BM25A). Among these 
constructs, BCoV DI RNA with a 1- and 2-nt deletion at the 3' 
terminus still supported the (+)-strand synthesis (54.9% and 84.7% 
of the wt BM25A, respectively) (Fig. 4C, lanes 8-9), whereas DI 
RNA constructs with a 3- and 4-nt deletion totally abolished the (+ 
)-strand DI RNA synthesis (Fig. 4C, lanes 6-7). Further deletion of 
5-, 15-, 30-, and 55- nts (constructs A5, A 15, A30 and A55, 
respectively) (Fig. 4A) from the 3' terminus of the 55-nt region led 
to a significant decrease in the synthesis of (-)-strand (less than 
10% of wt BM25A) (Fig. 4B) and (+)-strand DI RNA (Fig. 4C, 
lanes 2-5). 

To further dissect the specific cu-acting RNA motif required for 
(-)- and (+)-strand RNA synthesis within this 55-nt region, DI 
RNA constructs with various deletions from the 5' terminus of this 
region were created and tested (Fig. 4A). The efficiency of (-)- 
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Figure 2. Diagram of the hybrid BCoV Dl RNA and the head-to-tail strategy used to detect the synthesis of (-)-strand Dl RNA. (A) 

Upper panel: the composition of naturally occurring BCoV Dl RNA relative to the full-length BCoV genome. Lower panel: the modified BCoV Dl RNA 
construct BM 25A used in this study. The 65-nt leader sequence is illustrated with a filled rectangle. (B) Strategy for determining the synthesis of (-)- 
strand Dl RNA. RNA extracted at 8 hpt of Dl RNA transcript was treated with TAP and head-to-tail ligated. MHV-A59 3' UTR-negative-strand-specific 
primer 1 MHV3UTR3(-) was used for RT with the ligated RNA as template. MHV-A59 3' UTR-negative-strand-specific primer 2 MHV3UTR6(-) and BCoV 
5' UTR-positive-strand specific primer 3 BCoV23-40(+) were used for PCR. 
doi:1 0.1 371 /journal.pone.0098422.g002 



strand Dl RNA synthesis from Dl RNA with a deletion from nts — 
55 to -45 (construct A55-45) was not affected and was even better 
than that from wt BM25A (Fig. 4B) although the efficiency of (+)- 
strand Dl RNA synthesis was impaired (44.5% of wt BM25A, 
Fig. 4C, lane 10). The extension of the deletion from nts -55 to -40 
in Dl RNA (construct A55-40) still did not affect the efficiency of 
both (-)-and (+)-strand Dl RNA synthesis (Fig. 4 B and Fig. 4C, 
lane 1 1, respectively). Surprisingly, Dl RNA with a deletion from 
nts -55 to -35 (construct A55-35) almost blocked the synthesis of (— 
)-strand Dl RNA (less than 20% of wt BM25A, Fig. 4B) but still 
maintained the (+)-strand Dl RNA synthesis (78.9% of wt BM25A, 
Fig. 4C, lane 12). On the other hand, both (-)- and (+)-strand Dl 
RNA synthesis were significandy inhibited from Dl RNA with a 
deletion from nts -55 to -30 (construct A55-30) (Fig. 4B and 
Fig. 4C, lane 13, respectively). The efficiency of (-)-strand Dl RNA 
synthesis was also significantly blocked in the construct A30-15 
with a deletion from nts -30 to -15 (-40% of wt BM25A, Fig. 4B) 
and no (+)-strand Dl RNA synthesis was detected from this 
deletion mutant (Fig. 4C, lane 14), suggesting the nts from -30 to - 
15 are critical for both (-)- and (+)-strand Dl RNA synthesis. Of 
these deletion mutants, all of the mutated sequences in the 
synthesized (— )-strand Dl RNA were still maintained at 8 hpt; 
however, the available mutant sequences in the (-)- and (+)- strand 
Dl RNA at 48 hpi of VP1 showed that all contained the sequence 
with majority from wt Dl RNA except for the mutant A2 in which 
the sequence of mutant A2 was mixed with that of wt BM25A 
(Fig. 4C). Based on these results, the tentative conclusions at this 
time were as follows: (1) the nts from -1 to -4 and from -40 to -55 
are not important for the (-)-strand BCoV Dl RNA synthesis and 
(2) the nts from -1 to -2 and from -35 to -55 are not absolutely 
required for the (+)-strand BCoV Dl RNA synthesis. Therefore, 
the critical cu-acting elements essential for (-)- or (+)-strand 
synthesis may reside in the region between the nts -2 and -40. To 
further determine the requirement of cw-acting elements for the (— 



)- and (-t-)-strand BCoV Dl RNA synthesis in this region, deletion 
mutants A5-14, A15-24, A25-34, and A35-39 were constructed 
(Fig. 5A) and tested. As shown in Fig 5B, the efficiency of the (— )- 
strand Dl RNA synthesis from these deletion constructs was 
inhibited (-29% to -42% of wt BM25A). The (+)-strand Dl RNA 
synthesis was still maintained from construct A35-39 (60.4% of wt 
BM 25A) (Fig. 5C, lane 2) but was significantly blocked from 
constructs A5-14, A15-24 and A25-34 (Fig. 5C, lanes 3-5). 
Sequencing analysis showed that all of the mutated sequences in 
these synthesized (— )-strand Dl RNA were still maintained at 
8 hpt; however, all the mutants reversed back to wt Dl RNA at 
48 hpi of VP1. The sequencing results acquired from Dl RNA 
constructs Al, A55-45, A55-40, A55-35, and A35-39 at 48 hpi of 
VP1 (Figs. 4C and 5C) indicate that the (+)-strand Dl RNA 
detected by Northern blot assay may mostly originated from wt Dl 
RNA (BM25A) that likely arose by reversion and consequendy the 
importance of the mutated sequences remains to be determined. 
To confirm that these Dl RNA mutants are able to replicate, the 
supernatant collected at 48 hpt was used to infect fresh HRT-18 
cells. The total cellular RNA was harvested at 24 hpi of VP1 and 
then subjected to Northern blot assay and sequencing analyses. As 
shown in Fig. S2, the relative replication efficiency of these 
mutants at 24 hpi of VP1 was similar to that at 48 hpi of VP1. 
Sequencing analysis revealed that the mutated sequences were still 
maintained in the mutants A2, A55-40 and A35-39. Although the 
sequencing results presented a mixed population in the mutants 
Al, A55-45 and A55-35, the majority of the sequences were from 
the mutants rather than from the wt Dl RNA (BM25A). These 
results suggest that these Dl RNA mutants still maintain the ability 
to replicate although the replication efficiency among these 
mutants varied. Taken together, based on the deletion mutagen- 
esis experiments and sequence analysis, we conclude that the nts 
from —5 to —39 within the 3 '-terminal 55 nts are the ay-acting 
elements responsible for the (— )-strand BCoV Dl RNA synthesis 
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Figure 3. Analysis of the requirement of 3 -terminal 55 nts for the synthesis of (— )-strand BCoV Dl RNA. (A) Diagram of the BCoV Dl 
RNA BM25A with the intact 3'-terminal 55 nts and the mutant construct A55 with the deletion of 3'-terminal 55 nts (denotes with dashes). (B) 
Detection of (-)-strand BCoV Dl RNA with head-to-tail ligation and RT-PCR. RT-PCR products with a size of —150 bp were observed from BCoV- 
infected BM25A-transfected cells (lanes 2-8, arrowhead) but not from BCoV-infected A55-transfected cells (lanes 10-16). Lanes 18-21 represent the 
controls for RT-PCR. CI : total cellular RNA from mock-infected cells. C2: total cellular RNA from BCoV-infected cells. C3: total cellular RNA from Dl RNA- 
transfected mock-infected cells. C4: a mixture of BCoV-infected cellular RNA extracted at 10 hpi and 200 ng of BM25A transcript. (C) Sequence of the 
cDNA-cloned RT-PCR product with a size of —150 bp from lane 5 in Fig. 3B. [shown in the (+)-strand], (D) Detection of the potential recombination 
between the helper virus genome and Dl RNA. The primers MHV3' UTR2(+), which anneals to the MHV 3' UTR and was used for RT, and M3(-), which 
anneals to BCoV M protein gene, were used for PCR to detect potential recombination between helper virus BCoV genome and BM25A (lane 2) or 
A55 (lane 3). A recombinant DNA of 1,639 nt shown in lane 4 was created by overlap RT-PCR and was used as a size marker for the product generated 
with the primers MHV 3' UTR2(+) and M3(-). (E) Left panel: the relative levels of (-)-strand Dl RNA synthesis as measured by q RT-PCR. Control A: total 
cellular RNA from mock-infected cells. Control B: total cellular RNA from BCoV-infected cells. Control C: total cellular RNA from Dl RNA-transfected 
mock-infected cells. Control D: a mixture of BCoV-infected cellular RNA extracted at 1 0 hpi and 200 ng of BM25A transcript. Right panel: the amounts 
of Dl RNA, helper virus N sgmRNA, and 18S rRNA from Dl RNA-transfected BCoV-infected cells at 8 hpt of VP0 as measured by Northern blot assay. 
The values (E) represent the mean ±SD of three individual experiments. SD: standard deviation. ***p<0.001. 
doi:1 0.1 371 /journal.pone.0098422.g003 



and that the nts from —3 to — 34 within the 3'-terminal 55 nts are 
the ay-acting elements required for the (+)-strand BCoV Dl RNA 
synthesis. 



The effect of the 3'-most nucleotide species on the (-)- 
and (+)-strand Dl RNA synthesis 

Little is known regarding the mechanism of the initiation of (— )- 
strand RNA synthesis in coronaviruses. Since the last nucleotide of 
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Figure 4. Identification of as-acting RNA elements within the 3 -terminal 55 nts that are required for (— )- and (+)-strand RNA 
synthesis. (A) Constructs of deletion mutants within the 3'-terminal 55 nucleotides of BCoV Dl RNA. Dashes denote deleted sequences. (B) The 
relative levels of (-)-strand Dl RNA synthesis. Total cellular RNA was extracted from Dl RNA-transfected BCoV-infected cells at 8 hpt. The synthesis of 
(-)-strand Dl RNA from the deletion mutant was quantitated by qRT-PCR and was compared with that from wt BM25. Control A: total cellular RNA 
from mock-infected cells. Control B: total cellular RNA from BCoV-infected cells. Control C: total cellular RNA from Dl RNA-transfected mock-infected 
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cells. Control D: a mixture of BCoV-infected cellular RNA extracted at 10 hpi and 200 ng of BM25A transcript. (C) Upper panel: the synthesis of (+)- 
strand Dl RNA as detected by Northern blot assay. Total cellular RNA was extracted at 48 hpi of VP1 and was subjected to Northern blot assay with N 
sgmRNA and 1 8S rRNA as internal controls. Middle panel: the relative levels of the (+)-strand Dl RNA synthesis. Lower panel: the sequence of the BCoV 
Dl RNA at 48 hpi of VP1 as determined by RT-PCR and sequencing analysis. The values (B) and (C) represent the mean ±SD of three individual 
experiments. SD: standard deviation, wt: wild type, mx: mixed, NA: not available. *p<0.05, **p<0.01, ***p<0.001. 
doi:1 0.1 371 /journal.pone.0098422.g004 



the ?>' UTR in all coronavirus genomes sequenced to date is 
cytosine, it was speculated that this conserved nucleotide species 
may play an important role in coronavirus replication. To test this 
hypothesis, the cytosine residing at the last nucleotide of the 3' 
UTR in BCoV Dl RNA was substituted with adenine, uracil and 
guanine to make the constructs BM-3'A, BM-3'U, and BM-3'G, 
respectively (Fig. 6A). As shown in Fig. 6B, even though construct 
BM-3'A was still able to synthesize its (— )-strand counterpart 
(—60% of wt BM25A), its ability to synthesize (+)-strand was 
significantly abolished (7.5% of wt BM25A) (Fig. 6C, lane 2). In 
contrast, the efficiency of the (— )-strand Dl RNA synthesis from 
construct BM-3'U was significandy inhibited (~30% of wt 
BM25A) (Fig. 6B), but the efficiency of (+)-strand Dl RNA 
synthesis from construct BM-3'U was similar to that from wt 
BM25A (88.1% of wt BM25A, Fig. 6C, lane 3). The efficiency of 
both (-)- and (+)-strand Dl RNA synthesis (Fig. 6B and Fig. 6C, 
lane 4, respectively) from construct BM-3'G, on the other hand, 
was significantly blocked (-10% and -21.8% of wt BM25A, 
respectively). Sequencing analysis from VP1 revealed that the 



mutants BM-3'A and BM-3'G reverted back to the wt sequence, 
but the mutant BM-3'U still maintained the mutated sequence. 
The results together suggest that (1) the substitution of the 3 '-most 
nucleotide cytosine with other nucleotide species is sensitive for (— 
)-strand Dl RNA synthesis and (2) the substitution of the 3 '-most 
cytosine residue by uracil, but not adenine or guanine, is more 
tolerant for (+)-strand Dl RNA synthesis. The nucleotide species at 
the 3'-most position (— 1), therefore, is important, but not critical, 
for both (-)- and (+)-strand BCoV Dl RNA synthesis. 

Discussion 

Many as-acting elements required for coronavirus replication 
have been identified in the past two decades 
[3,4,5,6,7,8,9,10,11,12,13,14,15,16,17]. The identification of the 
Of-acting elements specific for (— )-strand RNA synthesis in 
coronavirus, however, has been hampered due to insufficiencies 
in the techniques used to detect and quantitate RNA (— ) strand 
species. In this study, we employed the method of head-to-tail 
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Figure 5. Determination of c/5-acting RNA elements between the nts —4 and —40 required for (— )- and (+)-strand RNA synthesis. (A) 

Deletion mutants of BCoV Dl RNA. Dashes denote deleted sequences. (B) The relative levels of (— )-strand Dl RNA synthesis. The synthesis of (— )- 
strand Dl RNA from the deletion mutant was quantitated by qRT-PCR and was compared with that from wt BM25. Control A: total cellular RNA from 
mock-infected cells. Control B: total cellular RNA from BCoV-infected cells. Control C: total cellular RNA from Dl RNA-transfected mock-infected cells. 
Control D: a mixture of BCoV-infected cellular RNA extracted at 10 hpi and 200 ng of BM25A transcript. (C) Upper panel: the synthesis of (+)-strand Dl 
RNA as detected by Northern blot assay. Total cellular RNA was extracted at 48 hpi of VP1 and was subjected to Northern blot assay with N sgmRNA 
and 18S rRNA as internal controls. Middle panel: the relative levels of the (+)-strand Dl RNA synthesis. Lower panel: the sequence of the BCoV Dl RNA 
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doi:10.1371/journal.pone.0098422.g005 
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ligation and qRT-PCR to examine the ay-acting elements required 
for (-)-strand BCoV Dl RNA synthesis [26,37]. This method has 
the advantages of detecting and quantitating low copy number of 
(— )-strand BCoV Dl RNA as well as distinguishing between T7 
RNA polymerase-generated (— )-strand Dl RNA transcripts and 
coronaviral polymerase-generated (— )-strand Dl RNA, and 
thereby allows us to systematically examine the ay-acting RNA 
elements required for (— )-strand Dl RNA synthesis and to gain 
insight into the mechanism of coronavirus replication. 

Replication of (+)-strand RNA viruses can be generally divided 
into two steps. First, the viral polymerase employs a (+)-strand 
genome as a template to generate the (— )-strand complement. 
Second, the synthesized (— )-strand genome then serves as a 
template for the synthesis of the (+)-strand counterpart. The 
majority of previous studies on ay-acting elements in coronaviruses 
emphasize their roles in replication [interpreted as (+)-strand 
synthesis]; however, it is still unclear that the identified ay-acting 
elements function in (— )- or (+)-strand synthesis. We in this study, 
therefore, simultaneously analyzed the effects of ay-acting elements 
on the synthesis of both (— )- and (+)-strand BCoV Dl RNA during 
infection in order to specifically define the role of a ay-acting 
element. According to the results concerning the (— )- and (+)- 
strand Dl RNA synthesis from mutants A5, A15, A30 and A55, we 



conclude that the deleted sequences may function at least in the 
(— )-strand Dl RNA synthesis because the deleted sequences 
significandy blocked the (— )-strand Dl RNA synthesis from their 
(+)-strand Dl RNA templates (Fig. 4B). Regarding the role of these 
deleted sequences in the (+)-strand Dl RNA synthesis, however, it 
is not feasible to speculate from these mutants at this time because 
the undetectable (+)-strand Dl RNA in VP 1 from these mutants 
may be due to the very low amounts of the (— )-strand Dl RNA 
templates or due to the deleted sequences which block the 
synthesis of (+)-strand Dl RNA. For the mutants Al, A2, A3, and 
A4, the (— )-strand Dl RNA synthesis was overall slightly affected 
with the deleted sequences (Fig. 4B). However, the defect in (+)- 
strand Dl RNA synthesis from the mutants A3 and A4 (Fig. 4C, 
lanes 6-7) indicates that nts from —3 to —4 are critical ay-acting 
elements in the synthesis of (+)-strand rather than (— )-strand Dl 
RNA. In comparison with wt BM25A, the efficiency of both (— )- 
and (+)-strand Dl RNA synthesis from the mutants A55-45 and 
A55-40 was not significandy inhibited, suggesting nts from —55 to 
—40 are not absolutely required for both (— )- and (+)-strand Dl 
RNA synthesis. Interestingly, the efficiency of (— )-strand RNA 
synthesis from deletion mutant A55-45 was even higher than that 
from wt BM25A as shown in Fig. 4B. Regarding the 3'-terminal 
55 nts in the MHV 3' UTR, the sequence between -1 and -44 
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in the deletion mutant A55-45 is nearly identical to that in the 
BCoV 3' UTR after alignment. It could be possible that BGoV 
RNA-dependent RNA polymerase utilizes the almost authentic 
BCoV 3' UTR in mutant A55-45 more efficiently than the MHV 
3' UTR in the BM25A, leading to a better efficiency in the 
synthesis of (— )-strand DI RNA. It is also intriguing that the 
replication efficiency from the 16-nts deletion mutant A55-40 was 
enhanced, but the replication efficiency from the 1 1 -nts deletion 
mutant A55-45 was repressed in comparison with that from wt 
BM25A (Fig. 4G). However, the involvement of 3' UTR elements 
such as enhancers and repressors in the regulatory functions of 
replication has been suggested in RNA viruses [38,39,40]. 
Therefore, considering the replication efficiency of A55— 45, 
A55-40 and A55-35 (44%, 102% and 78%, respectively), we 
speculate that nts —45 to —55 may serve as an enhancer because 
the deletion of these nts in mutant A55-45 leads to the modest 
inhibition of replication. In contrast, nts —40 to —44 may act as a 
repressor and the synergistic effect of both enhancer and repressor, 
which are missing in the mutant A55-40, could be the reason why 
the construct had the similar or even higher replication efficiency 
in comparison with wt BM25A. Similarly, nts —35 to —39 may 
play a role of enhancer as evidenced by the comparison of 
replication efficiency between A55-40 and A55-35 (102% vs 78%) 
(Fig. 4C) and between wt BM25A and A35-39 (100% vs 60%) 
(Fig. 5G). In mutants A55-35 and A35-39, the efficiency of (-)- 
strand DI RNA synthesis was significantly decreased, but the (+)- 
strand DI RNA synthesis was still maintained. The results suggest 
that the efficiency of (+)-strand DI RNA synthesis may not be 
affected by the existing numbers of (— )-strand DI RNAs as long as 
the cw-acting element [on (— )-strand DI RNA] required for (+)- 
strand DI RNA synthesis is intact. These results are also in 
agreement with the previously proposed model which indicated 
that the (— )-strand viral RNA is repeatedly used as a template for 
(+)-strand viral RNA synthesis [29]. Based on the criteria, the 
deleted sequences in mutants A5-14, A15-24 and A25-34 are 
required for both (— )- and (+)-strand DI RNA synthesis since the 
(— )- and (+)-strand DI RNA synthesis from these mutants was 
significantly inhibited (Figs. 5B and 5C). Therefore, simultaneously 
analyzing the efficiency of (— )- and (+)-strand DI RNA synthesis 
enables us to specifically define the Of-acting elements that 
function in (— )- or (+)-strand DI RNA synthesis. 

Because the last nucleotide of the 3' UTR in all coronavirus 
genomes sequenced to date is cytosine and coronavirus nonstruc- 
tural protein 8 has been proposed to preferentially initiate the 
synthesis of primer at the cytosine within the internal 5'-(G/ 
U)CC-3' trinucleotides [41], we in the present study speculated 
that the cytosine residue may play a critical role in (— )-strand 
RNA synthesis. Since the substitution of the cytosine with adenine, 
guanine or uracil inhibited the (— )-strand DI RNA synthesis at 
different levels (Fig. 6B), the 3'-most nucleotide species apparendy 
is critical for (— )-strand DI RNA synthesis. Nevertheless, whether 
the cytosine is involved in the primer-dependent (— )-strand RNA 
synthesis remains to be elucidated. On the other hand, the results 
that the substitution of the cytosine with adenine and guanine, but 
not uracil, almost blocked (+)-strand DI RNA synthesis suggest 
that a purine (adenine and guanine), but not a pyrimidine (cytosine 
and uracil), on the 5'-most nucleotide of (— )-strand 3' UTR is 
preferred for subsequent (+)-strand synthesis. Thus, these results 
suggest that the 3'-most nucleotide species on the DI RNA (+) 
strand is an important element for both (— )-and (+)-strand DI 
RNA synthesis. 

The initiation of (— )-strand RNA synthesis in (+)-strand RNA 
viruses occurs at the 3' terminal sequences of (+)-strand RNA, and 
consequendy alterations of the 3 '-terminal sequences of (-t-)-strand 



RNA may be detrimental to (— )-strand RNA synthesis as 
evidenced by this study and others [42,43,44,45,46,47,48]. 
However, evidence from this study and previous findings 
[12,13,14,15] has also shown that alterations of the sequence at 
the 3' terminal region of (+)-strand RNA [5' terminal region on (— 
)-strand RNA] also impairs the synthesis of (+)-strand RNA. The 
reasons why the sequence alterations at the 3 '-terminal region on 
(+)-strand RNA, for example, mutants A3, A4, BM-3'A, and EM- 
S' G, affected the synthesis of (+)-strand RNA remains unknown. 
However, since (1) the interactions between the 5' and 3' end of 
the genome facilitated by RNA or/and protein to form the 
replication complex and thereby to initiate replication have been 
well documented in numerous of RNA viruses [49,50,51,52,53,54] 
and (2) the results from this study that the mutations in the 3'- 
terminal region of (+)-strand DI RNA were able to alter the 
efficiency of (+)-strand DI RNA synthesis, we postulate that the 
sequence mutations in the 3 '-terminal region of (+)-strand DI 
RNA [5' terminal region of (— )-strand DI RNA] may hinder the 
binding of proteins to this region and then block the interactions 
between the 5' and 3' end of the (— )-strand genome to constitute 
the replication complex that is required for (+)-strand synthesis. 
This notion is supported by the previous studies in which the nts 
from —26 to —36 in the MHV 3' UTR have been demonstrated 
to be required for host protein binding and viral RNA replication 
[18,55]. Further study is required to elucidate whether the as- 
acting elements identified in this study and others are involved in 
the constitution of replication complex in the (— )-strand RNA 
genome, which is required for the subsequent (+)-strand RNA 
synthesis. 

In this study we performed the head-to-tail ligation of RNA and 
qRT-PCR to detect and quantitate coronavirus (— )-strand DI 
RNA synthesis. The method, along with Northern blot assay, 
allows us to simultaneously detect both (— )- and (+)-strand RNA 
synthesis and specifically define the functions of cu-acting elements 
in the (— )- or (+)-strand DI RNA synthesis. In coronaviruses, 
besides the as-acting elements identified here, many as-acting 
elements located in the 5' UTR, 3' UTR and coding region have 
been demonstrated to be required for replication; however, their 
specific function in (— )- or (+)-strand RNA synthesis remains to be 
determined. Using the method may contribute to elucidate the 
specific roles of these as-acting elements. 

Materials and Methods 

Viruses and cells 

A DI RNA-free stock of the Mebus strain of BCoV (GenBank 
accession no. U00735) at 3x10' PFU/ml was used as a helper 
virus in human rectal tumor (HRT)-18 cell line as described 
previously [56,57,58]. 

Plasmid constructs 

To generate the pBM25A construct, PGR was performed with 
the oligonucleotides TGEV 7(-) and BM25A(+) and with pDI 
RNA-1 in which the 288-nt 3' UTR of BCoV-Mebus in pDrepl 
was replaced with the 301-nt 3' UTR and 43-nt poly(A) tail of 
MHV-A59 (GenBank accession no. NC_001846) as a template 
[26,35]. The resultant PCR product was cloned into the TOPO- 
XL vector (Invitrogen) and digested with Spel and Mlul. The 
digested fragment was cloned into Spel- and MM-linearized pDI 
RNA-1 to create pBM25A containing MHV-A59 3' UTR and 25- 
nt poly(A) tail. 

To construct the mutant pA55, an overlap PCR mutagenesis 
procedure was performed as previously described [59] but with the 
oligonucleotides TGEV 7(-) and A55(+) and pBM25A DNA in 
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the first PCR, the oligonucleotides A55(— ) and DI reverse(+) and 
pBM25A DNA in the second PCR, and the oligonucleotides 
TGEV 7(— ) and DI reverse(+) and the products of the first two 
reactions in a third PCR to make a 1203-nt product that was 
cloned into the TOPO XL vector (Invitrogen). From this, a 756-nt 
fragment obtained by digestion with Spel and Mlul was cloned into 
Spel- and MM-linearized pBM25A to make the mutant pA55. 
Mutants of pAl, pA2, pA3, pA4, pA5, pAl5, pA30, pA55-45, 
pA55-40, pA55-35, pA55-30, pA30-15, pA5-14, pA 15-24, 
pA25-34, pA35-39, pBM-3'A, pBM-3'G, and pBM-3'U were 
similarly constructed, except for the corresponding oligonucleo- 
tides used in the first and second reactions as described in Table 
SI. 

Head-to-tail ligation and RT-PCR for detecting (-)-strand 
DI RNA products 

To synthesize transcripts in vitro, the MM-linearized plasmid 
DNA construct was transcribed with the RiboMAX Large Scale 
RNA Production System-T7 (Promega) according to the manu- 
facturer's instructions and was chromatographed through a 
Biospin 6 column (Bio-Rad) before use. To analyze the synthesis 
of the (-)-strand DI RNA (Fig. 2B), the HRT-18 cells in 35-mm 
dishes at ~80% confluency (~8xl0 5 cells/dish) were infected 
with BCoV at a multiplicity of infection of 5 PFU per cell. After 

2 hours of infection, 3 ug of BM25A or mutant transcript 
quantitated by RNA denaturing gel was transfected into BCoV- 
infected HRT-18 cells. The total cellular RNA was extracted with 
TRIzol (Invitrogen) at 8 hpt, and 10 ug of extracted RNA in 25 ul 
of water, 3 ul of 10X buffer and 10 U of (in 1 ul) tobacco acid 
pyrophosphatase (TAP) (Epicentre) were used to remove the 5'- 
triphosphate groups of the (-)-strand DI RNA, generating 
monophosphate DI RNA (-) strand. Following the treatment of 
TAP, phenol-chloroform-extracted RNA in 25 ul of water was 
heat-denatured at 95°C for 5 min and then quickly cooled. Then, 

3 ul of 10X ligase buffer and 2 U (in 2 ul) of T4 RNA ligase I 
(New England Biolabs) were added, and the mix was incubated for 
16 h at 16°C. Phenol-chloroform-extracted ligated RNA was 
quantitated, and 1 Ug of ligated RNA was used for RT reaction 
with oligonucleotide MHV3'UTR3(-), which binds nts 89-112 
from the poly(U) tail in the (-) strand of the MHV-A59 3' UTR, to 
synthesize cDNA with Superscript III reverse transcriptase 
(Invitrogen). Of this, 5 ul was used in a 50-ul PCR with 
AccuPrime Tag DNA polymerase (Invitrogen) and oligonucleotide 
MHV3'UTR6(-), which binds nts 64-85 from the poly(U) tail in 
the (-) strand of the MHV-A59 3' UTR, and BCV23-40(+), which 
binds nts 23—40 in the (+)-strand leader sequence of BCoV. The 
mixture was heated to 94°C for 2 min and then subjected to 50 
cycles of 30 s at 94°C, 30 s at 60°C, and 30 s at 72°C. The -150- 
bp PCR product was sequenced after cloning into TOPO XL 
vector (Invitrogen). 

Quantitation analysis of (-)-strand DI RNA synthesis by 
qRT-PCR 

For quantitating (-)-strand DI RNA synthesis, 1 ug of TAP- 
treated and ligated RNA was used for RT reaction with 
oligonucleotide MHV3'UTR3(-) and Superscript III reverse 
transcriptase (Invitrogen). Real-time PCR amplification using the 
primers MHV3'UTR6(-) and BCV23-40(+) was performed using 
TagMan Universal PCR Master Mix (Applied Biosystems) 
according to the manufacturer recommendations with a Light- 
Cycler 480 instrument (Roche Applied Science). To quantitate the 
synthesis of (— )-strand DI RNA, dilutions of plasmids containing 
the same sequence as the detected RT-PCR product of (-) -strand 



DI RNA were always run in parallel with quantitated cDNA for 
use as standard curves (dilutions ranged from 10 8 to 10 copies of 
each plasmid). In addition, 18S rRNA with primers 18S rRNA(-) 
and 18S rRNA(+), DI RNA with primers MHV3UTR2(-) and 
BCV29-54(+) and N sgmRNA with primers leader20(-) and 
BC VN(+) were applied as internal controls for normalization of (— 
)-strand DI RNA synthesis. The reactions were conducted with an 
initial pre-incubation at 95°C for 5 min, followed by 35 
amplification cycles as follows: 95°C for 15 s and 60°C for 60 s. 

Northern blot assay for DI RNA replication 

Ten micrograms of TRIzol-extracted total cellular RNA at 
48 hpi of VP1 was used for electrophoresis in a formaldehyde- 
agarose gel. RNA was transferred from the gel to Nytran 
membrane by vacuum blotting, and the blots were probed with 
oligonucleotide TGEV 8(+) (for DI RNA), BCVN(+) (for N 
sgmRNA) or 18SrRNA(+) (for 18S rRNA), which was tailed with 
digoxigenin (DIG)-ddUTP using a DIG Oligonucleotide 3 '-End 
Labeling kit (Roche Molecular Biochemicals). The RNA detected 
was visualized according to the manufacturer's recommended 
procedure. 

Sequence analysis of mutant DI RNA 

To identify the sequence of the (-)-and (+)-strand DI RNA at 
8 hpt and at 48 hpi of VP1, 10 ug of total cellular RNA was 
treated with tobacco acid pyrophosphatase (Epicentre) and ligated 
with T4 RNA ligase I (New England Biolabs). One microgram of 
ligated RNA was used for an RT reaction with oligonucleotide 
MHV 3' UTR3(-) [for (-)-strand DI RNA] or BCV107(+) [for (+)- 
strand DI RNA] to synthesize cDNA with Superscript III reverse 
transcriptase (Invitrogen), after which the oligonucleotides MHV 
3' UTR2(-) and BCV29-54(+) were used for PCR with AccuPrime 
Tag DNA polymerase (Invitrogen). The resultant PCR product 
was directiy sequenced. 

Supporting Information 

Figure SI Schematic diagram of the procedure to detect 
the synthesis of both (-)- and (+) -strand DI RNA. The 

HRT-18 cells were infected with helper virus BCoV and then 
transfected with DI RNA. The virus contained within transfected 
cells was referred to as virus passage 0 (VP0). The total 
intracellular RNA was collected at 8 hpt to detect the synthesis 
of (-)-strand DI RNA with qRT-PCR. The supernatant collected 
at 48 hpt from the BCoV-infected DI RNA-transfected cells 
contained virus passage 1 (VPl) and was used to infect fresh HRT- 
18 cells. The total intracellular RNA collected at 48 hpi of VPl 
was used to detect the synthesis of (+)-strand DI RNA with 
Northern blot assay. The 30-nt TGEV reporter is the binding site 
for the probe TGEV8(+) used in Northern blot assay. 
(TIF) 

Figure S2 The synthesis of (+)-strand DI RNA as 
detected by Northern blot assay. Total cellular RNA was 
extracted at 24 hpi of VPl and was analyzed by Northern blot 
assay with N sgmRNA and 18S rRNA used as internal controls. 
The sequence of the BCoV DI RNA at 24 hpi of VPl was 
determined by direct sequencing of RT-PCR products, wt: wild 
type, mut: mutant, mx: mixed. 
(TIF) 

Table SI Oligonucleotides used for this study. 

(TIF) 



PLOS ONE | www.plosone.org 



10 



May 2014 | Volume 9 | Issue 5 | e98422 



3'-Terminal Sequence for Negative- and Positive-Strand RNA Synthesis 



Acknowledgments 

We thank David A. Brian (University of Tennessee, Knoxville, TN) for 
many helpful discussions. 



References 

1. Woo PC, Lau SK, Huang Y, Yuen KY (2009) Coronavirus diversity, phylogeny 
and interspecies jumping. Exp Biol Med (Maywood) 234; 1117—1127. 

2. de Groot RJ, Baker SC, Barie R, Enjuancs L, Gorbalenya AE, et al. (2011) 
Family Coronaviridac. In: King AMGj Lefkowitz E, Adams MJ, Garstens EB, 
editors. Virus Taxonomy: Ninth Report of the International Committee on 
Taxonomy of Viruses: Elsevier Seicncc. pp. 806-828. 

3. Chang RY, Krishnan R, Brian DA (1996) The UCUAAAC promoter motif is 
not required for high-frequency leader recombination in bovine coronavirus 
defective interfering RNA. J Virol 70: 2720-2729. 

4. Kang H, Feng M, Schroeder ME, Giedroc DP, Leibowitz JL (2006) Putative cis- 
acting stem-loops in the 5' untranslated region of the severe acute respiratory 
syndrome coronavirus can substitute for their mouse hepatitis virus counterparts. 
J Virol 80: 10600-10614. 

5. Raman S, Bouma P, Williams CD, Brian DA (2003) Stem-loop III in the 5' 
untranslated region is a cis-acting element in bovine coronavirus defective 
interfering RNA replication. J Virol 77: 6720-6730. 

6. Raman S, Brian DA (2005) Stem-loop IV in the 5' untranslated region is a cis- 
acting element in bovine coronavirus defective interfering RNA replication. 
J Virol 79: 12434-12446. 

7. Li L, Kang H, Liu P, Makkinje N, Williamson ST, et al. (2008) Structural lability 
in stem-loop 1 drives a 5' UTR-3' UTR interaction in coronavirus replication. 
J Mol Biol 377: 790-803. 

8. Liu P, Li L, Keane SC, Yang D, Leibowitz JL, ct al. (2009) Mouse hepatitis virus 
stem-loop 2 adopts a uYNMG(U)a-like tetraloop structure that is highly 
functionally tolerant of base substitutions. J Virol 83: 12084-12093. 

9. Liu P, Li L, Millership JJ, Kang H, Leibowitz JL, ct al. (2007) A U-turn motif- 
containing stem-loop in the coronavirus 5' untranslated region plays a functional 
role in replication. RNA 13: 763-780. 

10. Gustin KM, Guan BJ, Dziduszko A, Brian DA (2009) Bovine coronavirus 
nonstructural protein 1 (p28) is an RNA binding protein that binds terminal 
genomic cis-replication elements. J Virol 83: 6087-6097. 

11. Brown CG, Nixon KS, Senanayake SD, Brian DA (2007) An RNA stem-loop 
within the bovine coronavirus nspl coding region is a cis-acting element in 
defective interfering RNA replication. J Virol 81: 7716-7724. 

12. Hsue B, Hartshorne T, Masters PS (2000) Characterization of an essential RNA 
secondary structure in the 3' untranslated region of the murine coronavirus 
genome. J Virol 74: 6911-6921. 

13. Hsue B, Masters PS (1997) A bulged stem-loop structure in the 3' untranslated 
region of the genome of the coronavirus mouse hepatitis virus is essential for 
replication. J Virol 71: 7567-7578. 

14. Liu Qj Johnson RF, Leibowitz JL (2001) Secondary structural elements within 
the 3' untranslated region of mouse hepatitis virus strain JHM genomic RNA. 
J Virol 75: 12105-12113. 

15. Williams GD, Chang RY, Brian DA (1999) A phylogcnetically conserved 
hairpin-type 3' untranslated region pseudoknot functions in coronavirus RNA 
replication. J Virol 73: 8349-8355. 

16. Guan BJ, Su YP, Wu HY, Brian DA (2012) Genetic evidence of a long-range 
RNA-RNA interaction between the genomic 5' untranslated region and the 
nonstructural protein 1 coding region in murine and bovine coronaviruses. 
J Virol 86: 4631-4643. 

17. Guan BJ, Wu HY, Brian DA (201 1) An optimal cis-replication stem-loop IV in 
the 5' untranslated region of the mouse coronavirus genome extends 16 
nucleotides into open reading frame 1. J Virol 85: 5593—5605. 

18. Yu W, LeibowitzJL (1995) A conserved motif at the 3' end of mouse hepatitis 
virus genomic RNA required for host protein binding and viral RNA replication. 
Virology 214: 128-138. 

19. Kim YN, Jcong YS, Makino S (1993) Analysis of cis-acting sequences essential 
for coronavirus defective interfering RNA replication. Virology 197: 53-63. 

20. Lin YJ, Lai MM (1993) Deletion mapping of a mouse hepatitis virus defective 
interfering RNA reveals the requirement of an internal and discontiguous 
sequence for replication. J Virol 67: 611 0-6 1 1 8. 

21. Goebel SJ, Miller TB, Bennett CJ, Bernard KA, Masters PS (2007) A 
hypcrvariablc region within the 3' cis-acting element of the murine coronavirus 
genome is nonessential for RNA synthesis but affects pathogenesis. J Virol 81: 
1274-1287. 

22. Spagnolo JF, Hogue BG (2000) Host protein interactions with the 3' end of 
bovine coronavirus RNA and the requirement of the poly(A) tail for coronavirus 
defective genome replication. J Virol 74: 5053-5065. 

23. Hofmann MA, Scthna PB, Brian DA (1990) Bovine coronavirus mRNA 
replication continues throughout persistent infection in cell culture. J Virol 64: 
4108-4114. 

24. Wu HY, Ozdarendeli A, Brian DA (2006) Bovine coronavirus 5 '-proximal 
genomic acceptor hotspot for discontinuous transcription is 65 nucleotides wide. 
J Virol 80: 2183-2193. 



Author Contributions 

Conceived and designed the experiments: WYL HYW. Performed the 
experiments: WYL TYK HYW. Analyzed the data: WYL HYW. 
Contributed reagcnts/materials/analysis tools: HYW. Wrote the paper: 
WYL HYW. 



25. Schcnborn ET, Micrendorf RC Jr (1985) A novel transcription property of SP6 
and T7 RNA polymerases: dependence on template structure. Nucleic Acids Res 
13: 6223-6236. ' 

26. Wu HY, Brian DA (2010) Subgenomic messenger RNA amplification in 
coronaviruses. Proc Natl Acad Sci USA 107: 12257-12262. 

27. Gorbalenya AE, Snijder EJ, Spaan WJ (2004) Severe acute respiratory syndrome 
coronavirus phylogeny: toward consensus. J Virol 78: 7863-7866. 

28. Snijder EJ, Bredenbeek PJ, Dobbe JC, Thiel V, ZiebuhrJ, et al. (2003) Unique 
and conserved features of genome and proteome of SARS-coronavirus, an early 
split-off from the coronavirus group 2 lineage. J Mol Biol 331: 991—1004. 

29. Lin YJ, Liao CL, Lai MM (1994) Identification of the cis-acting signal for minus- 
strand RNA synthesis of a murine coronavirus: implications for the role of 
minus-strand RNA in RNA replication and transcription. J Virol 68: 8131- 
8140. 

30. van der Most RG, Luytjes W, Rutjes S, Spaan WJ (1995) Translation but not the 
encoded sequence is essential for the efficient propagation of the defective 
interfering RNAs of the coronavirus mouse hepatitis virus. J Virol 69: 3744— 
3751. 

31. Fischer F, Stcgen CF, Koetzner CA, Masters PS (1997) Analysis of a 
recombinant mouse hepatitis virus expressing a foreign gene reveals a novel 
aspect of coronavirus transcription. J Virol 71:51 48-5 1 60. 

32. Baric RS, Fu K, Schaad MC, Stohlman SA (1990) Establishing a genetic 
recombination map for murine coronavirus strain A59 complementation groups. 
Virology 177: 646-656. 

33. Fu K, Baric RS (1994) Map locations of mouse hepatitis virus temperature- 
sensitive mutants: confirmation of variable rates of recombination. J Virol 68: 
7458-7466. 

34. Lai MM (1992) Genetic recombination in RNA viruses. Curr Top Microbiol 
Immunol 176: 21-32. 

35. Chang RY, Hofmann MA, Scthna PB, Brian DA (1994) A cis-acting function for 
the coronavirus leader in defective interfering RNA replication. J Virol 68: 
8223-8231. 

36. Wu HY, Ke TY, Liao WY, Chang NY (20 1 3) Regulation of Coronaviral Poly(A) 
Tail Length during Infection. PLoS One 8: c70548. 

37. Kc TY, Liao WY, Wu HY (2013) A Lcaderless Genome Identified during 
Persistent Bovine Coronavirus Infection Is Associated with Attenuation of Gene 
Expression. PLoS One 8: c82176. 

38. Zhang G, Zhang J, Simon AE (2004) Repression and derepression of minus- 
strand synthesis in a plus-strand RNA virus rcplicon. J Virol 78: 7619-7633. 

39. Zhang J, Simon AE (2005) Importance of sequence and structural elements 
within a viral replication repressor. Virology 333: 301-315. 

40. Zhang X, Kim CH, Sivakumaran K, Kao C (2003) Stable RNA structures can 
repress RNA synthesis in vitro by the brome mosaic virus replicase. RNA 9: 
555-565. 

41. Imbert I, Guillemot JC, Bourhis JM, Bussetta C, Coutard B, et al. (2006) A 
second, non-canonical RNA-dependent RNA polymerase in SARS coronavirus. 
EMBOJ 25: 4933-4942. 

42. Deiman BA, Koenen AK, Verlaan PW, Pleij CW (1998) Minimal template 
requirements for initiation of minus-strand synthesis in vitro by the RNA- 
dependent RNA polymerase of turnip yellow mosaic virus. J Virol 72: 3965- 
3972. 

43. Jacobson SJ, Konings DA, Sarnow P (1993) Biochemical and genetic evidence 
for a pseudoknot structure at the 3' terminus of the poliovirus RNA genome and 
its role in viral RNA amplification. J Virol 67: 2961-2971. 

44. Kao CC, SunJH (1996) Initiation of minus-strand RNA synthesis by the brome 
mosaicvirus RNA-dependent RNA polymerase: use of oligoribonucleotide 
primers. J Virol 70: 6826-6830. 

45. Miller WA, BujarskiJJ, Dreher TW, Hall TC (1986) Minus-strand initiation by 
brome mosaic virus replicase within the 3' tRNA-like structure of native and 
modified RNA templates. J Mol Biol 187: 537-546. 

46. Singh RN, Dreher TW (1997) Turnip yellow mosaic virus RNA-dependent 
RNA polymerase: initiation of minus strand synthesis in vitro. Virology 233: 
430-439. 

47. You S, Padmanabhan R (1999) A novel in vitro replication system for Dengue 
virus. Initiation of RNA synthesis at the 3'-end of exogenous viral RNA 
templates requires 5'- and 3 '-terminal complementary sequence motifs of the 
viral RNA. J Biol Chem 274: 33714-33722. 

48. Hardy RW, Rice CM (2005) Requirements at the 3' end of the sindbis virus 
genome for efficient synthesis of minus-strand RNA. J Virol 79: 4630—4639. 

49. Alvarez DE, Filomatori CV, Gamarnik AV (2008) Functional analysis of dengue 
virus cyclization sequences located at the 5' and 3'UTRs. Virology 375: 223- 
235. 



PLOS ONE | www.plosone.org 



11 



May 2014 | Volume 9 | Issue 5 | e98422 



3'-Terminal Sequence for Negative- and Positive-Strand RNA Synthesis 



50. Alvarez DE, Lodeiro MF, Luduena SJ, Pictrasanta LI, Gamarnik AV (2005) 
Long-range RNA- RNA interactions circularize the dengue virus genome. J Virol 
79: 6631-6643. 

51. Herold J, Andino R (2001) Poliovirus RNA replication requires genome 
circularization through a protein-protein bridge. Mol Cell 7: 581-591. 

52. Ooms M, Abbink TE, Pham C, Berkhout B (2007) Circularization of the HIV-1 
RNA genome. Nucleic Acids Res 35: 525.3-5261. 

5.3. Filomatori GV, Lodeiro MF, Alvarez DE, Samsa MM, Pietrasanta L, ct al. 
(2006) A 5' RNA element promotes dengue virus RNA synthesis on a circular 
genome. Genes Dev 20: 22.38-2249. 

54. Villordo SM, Gamarnik AV (2009) Genome cyclization as strategy for flavivirus 
RNA replication. Virus Res 139: 230-239. 

55. Yu W, LeibowitzJL (1995) Specific binding of host cellular proteins to multiple 
sites within the 3' end of mouse hepatitis virus genomic RNA. J Virol 69: 2016- 
2023. 



56. Lapps W, Hoguc BG, Brian DA (1987) Sequence analysis of the bovine 
coronavirus nucleocapsid and matrix protein genes. Virology 157: 47—57. 

57. King B, Brian DA (1982) Bovine coronavirus structural proteins. J Virol 42: 
700-707. 

58. Tompkins WA, Watrach AM, Schmale JD, Schultz RM, Harris JA (1974) 
Cultural and antigenic properties of newly established cell strains derived from 
adenocarcinomas of the human colon and rectum. J Nad Cancer Inst 52: 1101 — 
1110. 

59. Ozdarcndcli A, Ku S, Rochat S, Williams GD, Scnanayake SD, et al. (2001) 
Downstream sequences influence the choice between a naturally occurring 
noncanonical and closely positioned upstream canonical heptameric fusion motif 
during bovine coronavirus subgenomic mRNA synthesis. J Virol 75: 7362-7374. 



PLOS ONE | www.plosone.org 



12 



May 2014 | Volume 9 | Issue 5 | e98422 



